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Abstract 

We explore a nuclear physics resolution to the discrepancy between the predicted standard 
big-bang nucleosynthesis (BBN) abundance of 7 Li and its observational determination in 
metal-poor stars. The theoretical 7 Li abundance is 3 - 4 times greater than the observational 
values, assuming the baryon-to-photon ratio, r] wma , p , determined by WMAP. The 7 Li problem 
could be resolved within the standard BBN picture if additional destruction of A = 7 isotopes 
occurs due to new nuclear reaction channels or upward corrections to existing channels. This 
could be achieved via missed resonant nuclear reactions, which is the possibility we consider 
here. We find some potential candidate resonances which can solve the lithium problem and 
specify their required resonant energies and widths. For example, a 1~ or 2~ excited state of 
10 C sitting at approximately 15.0 MeV above its ground state with an effective width of order 
10 keV could resolve the 7 Li problem; the existence of this excited state needs experimental 
verification. Other examples using known states include 7 Be + 1 — > 10 B(18.80 MeV), and 
7 Be+d 9 B(16.71 MeV). For all of these states, a large channel radius (a > 10 fm) is needed 
to give sufficiently large widths. Experimental determination of these reaction strengths is 
needed to rule out or confirm these nuclear physics solutions to the lithium problem. 



1 Introduction 



Primordial nucleosynthesis continues to stand as our earliest probe of the universe based 
on Standard Model physics. Accurate estimates of the primordial abundances of the light 
elements D, 4 He and 7 Li within standard Big Bang Nucleosynthesis (BBN) [IHH] are crucial 
for making comparisons with observational determinations and ultimately testing the the- 
ory. Primordial abundances are also a probe of the early universe physics [6j. Currently, 
the theoretical estimates of D and He match the observational values within theoretical 
and observational uncertainties [21 [5] at the baryon-to-photon ratio determined by the 7- 
year WMAP data, r/ wmap = 6.19 ± 0.15 x 10~ 10 [TJ. In contrast, the theoretical primordial 
abundance of 7 Li does not match the observations. 

At ^ wm ap, the predicted BBN abundance of 7 Li ia3 [5] 



The observed 7 Li abundance is derived from observations of low-metallicity halo dwarf stars 
which show a plateau [8] in (elemental) lithium versus metallicity, with a small scatter 
consistent with observational uncertainties. An analysis [H] of field halo stars gives a plateau 
abundance of 



However, the lithium abundance in several globular clusters tends to be somewhat higher [TOj 
[TT] . and a recent result found in [11] gave 7 Li/H = (2.34±0.05) x 10~ 10 . Thus the theoretically 
estimated abundance of the isobar with mass 7 ( 7 Be+ 7 Li) is more than the observationally 
determined value by a factor of 2.2 - 4.2 [5], at f| wmap . Relative to the theoretical and 
observational uncertainties, this represents a deviation of 4.5-5.5 a. 

This significant discrepancy constitutes the "lithium problem" which could point to lim- 
itations in either the observations, our theoretical understanding of nucleosynthesis, or the 
post-BBN processing of lithium. 

On the theoretical front, strategies which have emerged to approach the lithium problem 
broadly either address astrophysics or microphysics. On the astrophysical side, one might 
attempt to improve our understanding of lithium depletion mechanisms operative in stellar 
models [12]. This remains an important goal but is not our focus here. 

The microphysical solutions to the lithium problem all in some way change the nuclear 
reactions for lithium production in order to reduce the primordial (or pre-Galactic) lithium 
abundance to observed levels. Some of these work within the Standard Model, focussing 
on nuclear physics, in particular the nuclear reactions involved in lithium production. One 
approach is to attempt to utilize the experimental uncertainties in the rates [2"] IT51 - H5" J. A 
second, related approach is the inclusion of new effects in the nuclear reaction database 
such as poorly understood resonance effects [16] . Finally, it may happen that effects beyond 

■"■Note that the 7 Li abundance reported here differ slightly from that given in [5], primarily due to the 
small shift in r\ as reported in (7]. 
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the Standard Model are responsible for the observed lithium abundance. For example, the 
primordial lithium abundance can be reduced by cosmological variation of the fine structure 
constant associated with a variation in the deuterium binding energy [17] , or by the post-BBN 
destruction of lithium through the late decays of a massive particle in the early universe [18J. 

In this paper, we remain within the Standard Model, examining the possible role of 
resonant reactions which may have been up to now neglected. The requisite reduction in the 
7 Li abundance can be achieved by either an enhancement in the rate of destruction of 7 Li 
or its mirror nucleus 7 Be. This approach is more promising than the alternative of reducing 
the production of 7 Be and 7 Li where the reactions are better understood experimentally and 
theoretically [T9H2T] . whereas the experimental and especially the theoretical situation for 
A = 8 — 11 has made large strides but still allows for surprises at the levels of interest to 
us [22]. 

The use of resonant channels is an approach that has paid off in the past in the context 
of stellar nucleosynthesis. Fred Hoyle famously predicted a resonant energy level at 7.68 
MeV in the 12 C compound nucleus which enhances the 8 Be + a — > 12 C reaction cross- section 
and allows the triple alpha reaction to proceed at relatively low densities [23]. Recently, 
it was shown that there are promising resonant destruction mechanisms which can achieve 
the desired reduction of the total A = 7 isotopic abundance [16] . This paper points to a 
resonant energy level at (E,J n ) = (16.71 MeV,5/2 + ) in the 9 B compound nucleus which 
can increase the rate of the 7 Be(d,p)aa and/or 7 Be(<i, 7) 9 B and thereby reduce the Be 
abundance. Here, we take a more general approach and systematically search for all possible 
compound nuclei [23] and potential resonant channels which may result in the destruction 
of 7 Be and/or 7 Li. 

Because of the large discrepancy between the observed and BBN abundance of 7 Li, any 
nuclear solution to the lithium problem will require a significant modification to the existing 
rates. As we discuss in the semi-analytical estimate in section [21 any new rate or modification 
to an existing one, must be 2 - 3 times greater than the current dominant destruction channels 
namely, 7 Li(p, a) a. for 7 Li and 7 Be(n,p) 7 Li for 7 Be. As discussed in [15] and as we show semi- 
analytically in § [21 this is difficult to achieve with non-resonant reactions. Hence, we will 
concentrate on possible resonant reactions as potential solutions to the lithium problem. As 
we will show, there are interesting candidate resonant channels which may resolve the 7 Li 
problem. For example, there is a possibility of destroying 7 Be through a l - or 2 _ 10 C excited 
state at approximately 15.0 MeV. The energy range between 6.5 and 16.5 MeV is currently 
very poorly mapped out and a state near the entrance energy for 7 Be + 3 He could provide 
a solution if the effective width is of order 10 keV. We will also see that these reactions all 
require fortuitously favorable nuclear parameters, in the form of large channel radii, as also 
found by Cyburt and Pospelov [16] in the case of 7 Be + d. Even so, in the face of the more 
radical alternative of new fundamental particle physics, these more conventional solutions 
to the lithium problem beckon for experimental testing. 

The paper is organized as follows: First, we lay down the required range of properties 
of any resonance to solve the lithium problem by means of a semi-analytic estimate in- 
spired by [23|[2^] in § [21 Then, in § [31 we list experimentally identified resonances from 
the databases: TUNL [24J and NNDC [21], involving either the destruction of 7 Be or 7 Li. 
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Finally, the solution space of resonant properties, wherein the lithium problem is partially 
or completely solved, is mapped for the most promising initial states involving either 7 Li or 
7 Be, by including these rates in a numerical estimation of the 7 Li abundance. This exercise 
will delineate the effectiveness of experimentally studied or identified resonances as well as 
requirements of possible missed resonant energy levels in compound nuclei formed by these 
initial states. This is described in § [4j We note that in our analysis, the narrow resonance 
approximation is assumed which may not hold true in certain regions of this solution space. 
Our key results are pared down to a few resonant reactions described in § |5j A summary 
and conclusions are given in § El 



2 Semi-analytical estimate of important reaction rates 

Before we embark on a systemic survey of possible resonant enhancements of the destruction 
of A = 7 isotopes, it will be useful to estimate the degree to which the destruction rates must 
change in order to have an impact on the final 7 Li abundance. The net rate of production 
of a nuclide % is given by the difference between the production from nuclides k and I and 
the destruction rates via nuclide j, i.e. for the reaction i + j — > k + l. This is expressed 
quantitatively by the rate equation [28] for abundance change 

-J- = -3Hm + y^n k ni{av) k i - UiUjiav)^, (3) 
jkl 

where is the number density of nuclide i, H is the Hubble parameter, riinj{av)ij are the 
sum of contributions from all the forward reactions destroying nuclide i and Ylki n k n i( av )ki 
are the reverse reactions producing it. (cry) is the thermally averaged cross-section of the 
reaction. The dilution of the density of these nuclides due to the expansion of the universe 
can be removed by re-expressing eq. ([3]) in terms of number densities relative to the baryon 
density Yi = rii/n h , as, 

dY 

— - = nb^^HH-r^Hy • (4) 
jkl 

Using this general form, the net rate of 7 Be production can be approximated in terms of the 
thermally averaged cross-sections of its most important production and destruction channels 
as 



dt 



rib ((crv)3 Uea Y3 He Y a - (av)7 Ben Y7 Be Y n ) . (5) 



Here, the reverse reaction rates of these production and destruction channels are neglected, 
as they are much smaller than the forward rates at the lithium synthesis temperature. A 
similar equation can be written down for 7 Li. When quasi-static equilibrium is reached, the 
destruction and production rates are equal. In this case, approximate values for new rates, 
which can effectively destroy either isobar, can be obtained analytically. 
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At temperatures T rs 0.04 MeV, both 7 Li and 7 Be are in equilibrium [26J which gives, 

(av) SRea Y 3Uc Y a = (av)r Bea Yr Be Y n . (6) 

Consider a new, inelastic 7 Be destruction channel 7 Be + X —*Y + Z, involving projectile X. 
This reaction will add to the right hand side of eq. fl6]) and shift the equilibrium abundance 
of 7 Be to a new value as follows, 

7Be ~ Mr Ben y n + (av)r BeX Y x ~ j + wyx 7 Be • u 

If the new reaction is to be important in solving the lithium problem, it must reduce the 
7 Be abundance by a factor of Yr^/Yf^ ~ 3 — 4 . This in turn demands via eq. (JTj) that 
(crv)7 BcX Yx/ (crv)7 Bcn Y n ~ 2 — 3, i.e., the rate for the new reaction exceeds that of the usual 
n — p interconversion reaction rate. A similar estimate can be made for 7 Li. 

This reasoning would exclude non-resonant rates as they would be required to have 
unphysically large astrophysical S-factors in the range of order 10 5 — 10 9 keV - barn depending 
on the channel. Thus we would expect that only resonant reactions can produce the requisite 
high rates. Possible resonant reactions are listed in the next section, whose key properties of 
resonance strength, r e g and energy, E res , lie in appropriate ranges capable of achieving the 
required destruction of mass 7. 

Finally, we turn to 7 Li destruction reactions, 7 Li+X — > Y + Z. Recall that at the WMAP 
value of r), mass 7 is made predominantly as 7 Be, with direct 7 Li production about an order 
of magnitude smaller. This suggests that enhancing direct 7 Li destruction will only modestly 
affect the final mass-7 abundance; we will see that this expectation is largely correct^] 

With these pointers, the list in the next section is reduced and numerical analysis of the 
remaining promising rates is done. 



3 Systematic Search for Resonances 

In this section we describe a systematic search for nuclear resonances which could affect 
primordial lithium production. We first begin with general considerations, then catalog the 
candidate resonances. We briefly review the basic physics of resonant reactions to establish 
notation and highlight the key physical ingredients. 

2 A subtle point is that normally, the mass-7 abundance is most sensitive to rate 7 Be(n,p) 7 Li [25] , 
Of course, this reaction leaves the mass-7 abundance unchanged, but the lower Coulomb barrier for 7 Li 
leaves it vulnerable to the 7 Li(p, a) 4 He reaction, which is extremely effective in removing 7 Li. Thus, for a 
new, resonant 7 Li destruction reaction to be important, it must successfully compete with the very large 
7 Li(p, a) 4 He rate, and even then the mass-7 destruction "bottleneck" remains the 7 Be(n,p) 7 Li rate that 
limits 7 Li appearance. Thus we would not expect direct 7 Li destruction to be effective. We will examine 7 Li 
destruction below, and confirm these expectations. 
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3.1 General Considerations 



Energetically, the net process 7 Be + A — > B + D must have Q + E- m [ t > 0, where the initial 
kinetic energy E^t ~ T < 40 keV is small at the epoch of A = 7 formation. Thus we in 
practice require exothermic reactions, Q > 0. Moreover, inelastic reactions with large Q will 
yield final state particles with large kinetic energies. Such final states thus have larger phase 
space than those with small Q and in that sense should be favored. 

Consider now a process 7 Be + X — > C* — > Y + Z which destroys 7 Be via a resonant 
compound state; a similar expression can be written for 7 Li destruction. In the entrance 
channel 7 Be + X — > C* the energy released in producing the compound state is Qc = 
A( 7 Be) + A(X) — A(C g ' s '), where A(A) = m — Am n is the mass defect. If an excited state 
C* in the compound nucleus lies at energy E ex , then the difference 

-Ercs = E cx — Q c (8) 

determines the effectiveness of the resonance. We can expect resonant production of C* 
if E vcs < T. In an ordinary ( "superthreshold" ) resonance we then have E ves > 0, while a 
subthreshold resonance has E res < 0. 

Once formed, the excited C* level can decay via some set of channels. The cross section 
for 7 Be + X — > C* — > Y + Z is given by the Breit-Wigner expression 

a[ ~ 2fiE (E — E rcs ) 2 + (r tot /2) 2 1 ' 

where E is the center-of-mass kinetic energy in the initial state, [i is the reduced mass and 

2J C , + 1 

W - (2J X + 1)(2J 7 +1) (10) 

is a statistical factor accounting for angular momentum. The width of the initial state 
(entrance channel) is T^it, and the width of the final state (exit channel) is Tfin. 

One decay channel which must always be available is the entrance channel itself. Obvi- 
ously such an elastic reaction is useless from our point of view. Rather, we are interested 
in inelastic reactions in which the initial 7 Be (or 7 Li) is transformed to something else. In 
some cases, an inelastic strong decay is possible where the final state particles Y + Z are 
both nuclei. Note that it is possible to produce a final-state nucleus in an excited state, e.g., 
C* — > Y* + Z, in which case the energy release Q' c is offset by the Y* excitation energy. 
This possibility increases the chances of finding energetically allowable final states. Indeed, 
such a possibility has been suggested in connection with the 7 Be + d — > 9 B* — > 8 Be* + p 
process [16J. 

Regardless of the availability of a strong inelastic channel, an electromagnetic transition 
C* — > (7^ +7 to a lower level is always possible. However, these often have small widths and 
thus a small branching ratio T^ n /T tot . Thus for electromagnetic decays to be important, a 
strong inelastic decay must not be available, and the rest of the reaction cross section needs 
to be large to compensate the small branching; as seen in eq. ()9]), this implies that r ini t be 
large. 
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Note that in all charged-particle reactions, the Coulomb barrier is crucially important 
and is implicitly encoded via the usual exponential Gamow factor in the reaction widths of 
both initial and final charged-particle states. However, if the reaction has a high Q, the final 
state kinetic energy will be large and thus there will not be significant final-state Coulomb 
supression; this again favors final states with large Q. In addition, if the entrance or exit 
channel has orbital angular momentum L > 0, there is additional exponential suppression, 
so that L > states are disfavored for our purposes. 

With these requirements in mind, we will systematically search for resonant reactions 
which could ameliorate or solve the lithium problem. We begin by identifying possible 
processes which are 

1. new resonances not yet included in the BBN code; 

2. 2-body to 2-body processes, since 3-body rates are generally very small in BBN due 
to phase space suppression as well as the relatively low particle densities and short 
timescales; 

3. experimentally allowed - in practice this means we seek unidentified states in poorly 
studied regimes; 

4. narrow resonances having r tot < T, which is around r tot < 40 keV but we will also 
consider somewhat larger widths to be conservatively generous. 

5. relatively low-lying resonances with E rcs < few x T ~ 100 — 300 keV, which are 
thermally accessible; here again we err on the side of a generous range. 

Once we have identified all possible candidate resonances, we will then assess their viability 
as solutions to the lithium problem based on available nuclear data. 

3.2 List of Candidate Resonances 

As described above, we will explore the resonant destruction channels of both 7 Li and 7 Be. 
Some of the potential resonances which might be able to reduce the mass 7 abundance to the 
observed value were recently considered in [16] . This analysis eliminates several candidate 
resonances, leaving as genuine solutions only the resonance related to the 7 Be(<i, 7) 9 B and 
7 Be(d, p)aa reactions and associated with the 16.71 MeV level in the 9 B compound nucleus. 
Here, we make an exhaustive list of possible promising resonances that may be important 
to either 7 Be or 7 Li destruction channels. In order to do so systematically and account for 
all possible resonances that may be of importance, we study the energy levels in all possible 
compound nuclei that may be formed in destroying 7 Be or 7 Li, making extensive use of 
databases at TUNL and the NNDC [SHETj. 

The available 2-body destruction channels 7 A+X may be classified by X = n, p, d, t, 3 He, a, 
and 7. Consequently, the compound nuclei that can be formed starting from mass 7 have 
mass numbers ranging from A = 8 to A = 11, and the ones of particular interest are 
8 Li, 8 Be, 8 B, 9 Be, 9 B, 10 Be, 10 B, 10 C, n B and n C. All relevant, resonant energy levels in these 
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compound nuclei that may provide paths for reduction of mass-7 abundance are listed in 
Tables HJ-HD1 



There are quantum mechanical and kinematic restrictions to our selection of candidates. 
The candidate resonant reactions must obey selection rules. The partial widths for a channel, 
which may be viewed as probability currents of emission of the particle in that channel 
through the nuclear surface, are given as 



where a is the channel radius and E is the projectile energy. Here k is the wavenumber of the 
colliding particles in the centre-of-mass frame and j 2 is the reduced width, which depends 
on the overlap between the wavefunctions inside and outside the nuclear surface, beyond 
which the nuclear forces are unimportant. The reduced width, 7 2 is independent of energy 
and has a statistical upper limit called the Wigner limit given by [30] 



The pre-factor of | is under the assumption that the nucleus is uniform and can change 
to within a factor of order unity if this assumption changes. The Wigner limit depends 
sensitively on the channel radius and thus varies with the nuclei involved. For the nuclei of 
our interest, typical values of 7 2 range from a few hundred keVs to a few Me Vs. 

In eq. (ITT)) . Pl(E,o) is the Coulomb penetration probability for angular momentum L 
and is a strong and somewhat complicated function of E and a. Thus, while the Wigner 
limit sets a theoretical limit on the reduced width, the upper limit on the full width, F L (E), 
depends on the values of Pl(E,o) and is sensitive to the details of the resonant channel 
being considered. In light of this complexity, our strategy is as follows. We evaluate the 
Tl(E) needed to make a substantial impact on the lithium problem. Then for the cases of 
highest interest, we will compare our results with the theoretical limit set by the Coulomb 
suppressed Wigner limit for those specific cases. 

We also limit our consideration to two body initial states, with resonance energies E res < 
650 keV. The high resonance energy limit ensures that all possible resonances which may 
influence the final 7 Li abundance are taken into account, though many of the channels with 
such high resonance energies will inevitably be eliminated. Excited final states have also been 
considered in making this list. Different excited states of final state products are marked as 
separate entries in the table, since each one has its own spin and therefore a different angular 
momentum barrier. And thereby the significance of each excited state in destroying mass 
7 is varied. Also, we usually eliminate the reactions with a negative Q-value except for the 



7 U{d,p) 8 U, 7 Be(d, 3 He) 6 Li, 7 Be(d,p) 8 Be* (16.922 MeV) and 7 Li( 3 He,p) 9 Be* (11.283 MeV) 



as they are only marginally endothermic. 

For a number of the reactions listed in these tables, 1-10, the total spin of the initial 
state reactants is equal to that of the compound nucleus, which is equal to the total spin of 
the products, with L = 0. However, for many reactions, angular momentum is required in 
the initial and/or final state, which decreases the penetration probability and thereby the 




(11) 
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width for that particular channel. In fact for some of these reactions, parity conservation 
demands higher angular momentum which worsens this effect. However, we do not reject any 
channel based on the angular momentum suppression of its width in these tables. Later we 
will shortlist those resonant channels which are not very suppressed and indeed potentially 
effective in destroying mass 7. 

The reactions of interest are listed in increasing order of the mass of the compound nuclei 
formed. The particular resonant energy levels of interest E ex and their spins are listed in the 
table. In general, different initial states involving 7 Li and 7 Be can form these energy levels 
and so all these relevant initial states are listed. For each one, the various final product 
states for an inelastic reaction are enumerated. Again, each of the final state products can 
also be formed in an excited state. These excited states must have lower energy than the 
initial state energies for the reaction to be exothermic. In addition spin and parity must 
be conserved. Enforcing these, the minimum final state angular momenta Lg n are evaluated 
from the spin of the resonant energy level and are listed in the tables. The total widths of 
the energy levels are listed whenever available. The partial widths of the different channels 
including the elastic one, out of each energy level are also listed. 

We adopt the narrow resonance approximation to evaluate the effect of these resonances 
and either retain or dismiss them as potential solutions to the lithium problem. Some of 
the partial widths or limits on them are high enough that they easily qualify to be broad 
resonances. This implies that the narrow resonance formula used to see their effect is not 
precise, but still gives a rough idea of whether the resonance is ineffective or not. 

Our expression for thermonuclear rates in the narrow resonance approximation is ex- 
plained in detail in Appendix |AJ and is given by 



This rate is controlled by two parameters specific to the compound nuclear state: E ves and 
r e ff. Here E res is given in eq. ©, and measures the offset from the entrance channel and the 
compound state. The resonance strength is quantified via 



with Tinit and Ta n being the entrance and exit widths of a particular reaction, and r to t the 
sum of the widths of all possible channels. Of these widths, the smaller of r init and Ta n along 
with r tot are listed in the table above. The resonance strength, r e ff ~ r init , if Tfi n dominates 
the total width, and vice versa. If r init and T^ n are the dominant partial widths and they 
are comparable to each other, then the strength is even higher. 

As discussed in Appendix |A] our narrow resonance rate in eq. f JTB"]) improves on the form 
of the usual expression for narrow resonance in two ways: (a) it extends to the subthresold 
domain; and (b) it introduces the factor / which accounts for a finite E res /T tot ratio. 

It is important to make a systematic and comprehensive search for all possible experimen- 
tally identified resonances capable of removing this discrepancy. In addition, it is possible 
that resonances and indeed energy levels themselves were missed, especially at the higher 




3/2 



(13) 




(14) 
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energies, where uncertainties are greater. Therefore it is useful to map the parameter space 
where the lithium discrepancy is removed to apriori lay down our expectations of such missed 
resonances. This can be done by looking at interesting initial states involving 7 Li and 7 Be, 
and abundant projectiles p, n, d, t, 3 He, a, and parametrizing the effect of inelastic channels 
on the mass- 7 abundance. This is described in § HI 

4 Narrow resonance solution space 

In order to study the effect of resonances in different compound nuclei on the abundance of 
mass 7, our strategy is to specify the reaction rate for possible resonances, and then run the 
BBN code to find the mass-7 abundance in the presence of these resonances. In particular, 
for reactions involving light projectile X, we are interested in considering the general effect 
of states 7 A + X — > C*, including those associated with known energy levels in the compound 
nucleus, as well as possible overlooked states. 

We assume that the narrow resonance approximation holds true at least as a rough guide. 
If the reaction pathway is specified, i.e., all of the nuclei 7 A+X — > C* — > Y + Z are identified, 
then the reduced mass fi, reverse ratio and the Q-value are uniquely determined. In this 
case, the thermally averaged cross-section is given by eq. (fl3l) . with two free parameters: the 
product uT e ff and the resonance energy, E res . Because the state C* is unspecified, so is its 
spin J*. On the other hand, we do know the spins of the initial state particles, and thus u 
is specified up to a factor 2J* + 1 (eq. [10]). For this reason, the Co»r eff dependence reduces to 
(2 J* + l)r cff , which we explicitly indicate in all of our plots. 

In a few cases we will be interested in one specific final quantum state, e.g., 7 Be(t, 3 He) 7 Li; 
when the final state is specified, the reaction can be completely determined, including the 
effect of the reverse rates. However, in most situations we are interested in the possibility 
of an overlooked excited state in the compound nucleus, and thus in unknown final states. 
In this scenario we thus have only a "generic" inelastic exit channel. Consequently, for such 
plots we cannot evaluate the reverse reaction rate (which is in all interesting cases small) 
and so we set the reverse ratio to zero. 

The resonant rates are included in the BBN code, individually for compound nuclei with 
an interesting initial state. The plots below show contours of constant, reduced mass-7 
abundances. A general feature of all the plots, is the near linear relation between logl^g 
and -E^es in the region of larger, positive values of T cS and E rcs . This can seen quantitatively 
as follows. The thermal rate is integrated over time or equivalently temperature to give the 
final abundance of mass 7 or 7 Li as it exists. Now assuming that the thermal rate operates at 
an effective temperature, Tu, at which 7 Li production peaks, a given value for this effective 
(av) will give a fixed abundance. This implies, 

5Y 7 /Y 7 ~ (av) peak ~ r eff e - £rcs/TLi ~ constant (15) 

This gives a feel for the linear relation in the plot. 



9 



Compound Nucleus, 


Initial 






-Eres 


Ttot 


Exit 


Exit Channel 




State 










( , n fi n ti pl 

\J lldllllClo 


Width 

v v iu lii 


8 Li, 3+, 2.255 MeV 


7 Li + n 


1 


1 


222.71 keV 


33 ± 6 keV 


7(ground state) 


7.0 ± 3.0 x 10~ 2 eV 


(Included) 






1 






n (elastic) «100% 


33 ± 6 keV 


8 Be, 2+ 16.922 MeV 


7 Li + p 


1 


2 

1 

2 
I 


-333.1 keV 


74.0 ± 0.4 keV 


7 (ground state) 
7(3.04 MeV) 
a « 100% 

"n 1 pl a p 1 

LJ l ClQG \jX\j J 


8.4 ± 1.4 x lO- 2 eV 
< 2.80 ±0.18 eV 
ss 74.0 keV 

U11IY11U W 11 


8d„ i + i 7 a,\c\ A/T^V 
r5e, i .D4U ivie v 


~Tt^ 

JjI + p 


1 
1 


i 
i 

1 

2 

1 

i 
i 


oo4.y Kev 


iu. * Ke v 


^(ground state) 
7(3.04 MeV) 
7(3.04 MeV) 
7(16.63 MeV) 
'vflfi 99 Mp~V) 

f \ XU.fti 1V1C V J 

p ^elastic j yo.o/o 


lo. ( ev 
6.7 ± 1.3 eV 
0.12 ±0.05 eV 
(3.2 ±0.3) x 10- 2 eV 

1 _L .U —1— U.OI A IU C V 

i u . o * Ke v 


8 Be, 2~ 18.91 MeV 


"'Be + n 





1 


10.3 keV 


122 keV* 


7(16.922 MeV) 


9.9 ± 4.3 x 10-' z eV 


(Included) 






1 

o 

Z 








7 (16.626 MeV) 
P 

p + -Li (0 iviev ) 
n (elastic) 


0.17 ±0.07 eV 

<r 1 05 1 keV* 
<. iuo.1 Kev 
16.65 keV* 


8 Re 3+ 1Q 07 MeV 


7 Be + 7i 


I 




1 70 3 keV 


970 4. on i,pV 

_ 1 \J I Zj w 1YC V 


n S3 1 00%, 


<r 970 keV 

Zi i VJ IvC V 


(Included) 






3 
1 
1 






p + 7 Li* (0.4776 MeV) 
7(3.03 MeV) 
n (elastic) 


< 270 keV 
10.5 eV 
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8 Re 3+ 19 935 MeV 
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I 


I 


335 3 keV 
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i 
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7(3.03 MeV) 
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10^ pV 
iu.o e v 

S3 1 1 3 5 keV 
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n 
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1 
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P 

p + 7 Li* (0.4776 MeV) 
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unknown 
unknown 
unknown 
unknown 


8 B S! ", 2+ MeV 


7 Be±p 


1 


1 


-0.1375 MeV 


unknown 


p (elastic) 
EC-> 8 Be 


unknown 















8.5 x 10- 19 eV 


8 B, 1+, 0.7695 MeV 


7 Be + p 


1 


1 


630 ±3 keV 


35.7 ±0.6 keV 


7 (ground state) 


25.2 ± 1.1 meV 


(Included) 






1 






p (elastic) 100% 


35.7 ±0.6 keV 



Table 1: This table lists the potential resonances in 8 Li, 8 Be and 8 B which may achieve 
required destruction of mass 7. These are all allowed by selection rules and includes some 
resonances already accounted for in determining the current theoretical 7 Li abundance in- 
dicated as (Included). The entrance and exit channels along with their partial and total 
widths (r tot ), minimum angular momenta (L init , L fin ) as well as resonance energies are listed 
wherever experimental data are available. The starred widths are a result of fits from R- 
matrix analysis. The list includes final products in ground and excited states with the latter 
marked with a star in the superscript. 

4.1 A = 8 Compound Nucleus 

As seen in Table [U the only resonance energy level of interest in the 8 Li compound nucleus 
at 2.255 MeV is already accounted for in the 7 Li + n reaction. In the 8 Be compound nucleus, 
there are six levels of relevance for destroying either 7 Li or 7 Be at 16.922, 17.64, 18.91, 19.07, 
19.24 and 19.40 MeV within our limit on E res . The 16.922 MeV level is more than 300 keV 
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below threshold and has a maximum total width of only 74 keV. Therefore, it is expected 
to have a weak effect. The 17.64 MeV level has typically low photon widths (~ 20 eV ) and 
a total width of 10.7 keV. But this state's decay is dominated by the elastic channel which 
makes this channel uninteresting. 

The energy level diagram for 8 Be [31] shows the initial state, 7 Be + n at an entrance 
energy of E = 18.8997 MeV bringing the 18.91, 19.07, 19.235 and 19.40 MeV levels into 
play. From among these the effect of the 18.91, 19.07 and 19.235 MeV resonances are 
already accounted for in the well known 7 Be(n,p) 7 Li reaction [19]. The (18.91 MeV, 2~), 
resonance with Lj n i t = is the dominant contributor [32~|l3"3] . Being a broad resonance with 
a total width of ~ 122 keV, the Breit-Wigner form is not used and instead an R- matrix fit to 
the data [19], is used to evaluate the contribution of the resonant rate. The remaining level 
at 19.40 should also contribute to this reaction through ground and excited states. Only 
the 19.40 MeV channel can have an a exit channel due to parity considerations. And this 
resonance, despite a high resonance energy of ~ 500 keV, can in principle be important due 
to its large total width of 645 keV, if the proton branching ratio is high. 

Figure [1] shows the 7 Li abundance in the (T e g , E res ) plane for the 7 Be(n,p) 7 Li reaction. 
Contours for 7 Li/H xlO 10 = 1.23, 2.0, 3.0, 4.0, and 5.0 (as labelled) are plotted as functions 
of the effective width and resonant energy. Below ~ (2J* + 1)40 = 120 keV, we expect 
our results based on the narrow resonance approximation to be quite accurate. As one can 
see from this figure, to bring the 7 Li abundance down close to observed values, one would 
require a very low resonance energy (of order ±30 keV) with a relatively large effective width. 
Unfortunately, the 19.40 MeV level of 8 Be corresponds to E Tes = 500 keV as shown by the 
vertical dashed line and does not make any real impact on the 7 Li abundance. 

Figure |2] shows the effect of a 8 B resonance with 7 Be and p in the initial state, plotted 
in the (r eff , -Eres) plane with contours of constant mass 7 abundances. According to Fig. [2] 
for resonance energies of a few 10 's of keV's, resonance strength of a few meV is sufficient 
to attain the observational value of mass 7. However, from the energy level diagram for 
8 B, [31], the closest resonant energy level, E ex is at 0.7695 MeV [34J, whose effect is already 
included via the 7 Be(p, 7) 8 B reaction. The experimental value of resonance energy is 632 
keV which is off the scale in this figure. The only other close energy level to the 7 Be + p 
entrance channel is at -0.1375 MeV which means that the ground state is a sub-threshold 
state. This is not the usual resonant reaction, since the ground state doesn't have a width in 
the sense we refer to a width for the other reactions. But at these energies, the astrophysical 
S"-factor is ~ 10 eV-barn which is very small and will yield a low cross-section. This too is 
off scale in the figure and verifies that the 7 Be(p, 7) 8 B reaction does not yield an important 
destruction channel. 

4.2 A = 9 Compound Nucleus 

The energy level diagram for 9 Be, [35J shows energy levels of interest at 16.671, 16.9752 and 
17.298 MeV; these appear in Table EJ The 7 Li + d entrance channel sits at 16.6959 MeV. The 
lowest lying resonant state is at 16.671 MeV and is a sub-threshold state with E res = —24.9 
keV which lies within the total width of 41 keV. This resonance is thus obviously tantalizing- 



11 




-100 100 200 300 400 500 600 
Resonance Energy E res [keV] 



Figure 1: The effect of resonances in the 8 Be compound nucleus involving initial states 
7 Be+n. It shows the range of values for the product of the resonant state spin degeneracy 
and resonance strength, (2 J + l)r e g, versus the resonance energy. Contours indicate where 
the lithium abundance is reduced to 7 Li/H = 1.23 x 10~ 10 (red), 2.0 x 1(T 10 (blue), 3.0 x 10~ 10 
(green) 4.0 x 10~ 10 (black) and 5.0 x 10~ 10 (magenta). Normal resonances have E res > 0, 
while subthreshold resonances lie in the E Tes < 0. The horizontal dot-dashed black line is the 
experimental value of the strength of the resonance corresponding to the 19.40 MeV energy 
level. The vertical dashed black line shows the position of E Tes for the same state. 
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Figure 2: As in Fig. [T]for the resonances in the 8 B compound nucleus involving initial states 
7 Be+p. 
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it is well-tuned energetically and involves an abundant, stable projectile. The 7 Li abundance 
contours for the 9 Be resonance states are shown in Fig. [31 Perhaps disappointingly, the 
figure shows that the effect on primordial mass 7 is minor. This illustrates the inability of 
direct 7 Li destruction channels to reduce the mass-7 abundance, as explained in §2j Given 
the overall difficulty of this channel, it is clear that the other possible resonant energy levels 
(16.9752 MeV and 17.298 MeV) also fail to substantially reduce the mass-7 abundance. 
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Table 2: As in Table [T] listing the potential resonances in 9 Be. 

The 9 B compound nucleus is relevant for studying the effect of the 7 ~Be(d,p)2a and its 
competitors such as 7 Be(d, 3 He) 6 Li and 7 Be(<i, a) 5 Li. As seen in Tabled the only two levels 
of interest here are the 16.71 and 17.076 MeV levels. The 16.71 MeV level corresponds to a 
resonance energy of 220 keV as shown by the vertical dashed line in Fig. H] The widths 
are unknown experimentally. The approximate narrow resonance limit on the resonance 
width which is shown by the horizontal solid line is around 40 keV. The p exit channel 
leads to the 7 Be(<i,p) 8 Be* reaction through the excited state at 16.63 MeV in 8 Be. This 
should eventually lead to formation of alpha particles. Fig. H] shows the effect of the 16.71 
MeV resonance on the mass-7 abundance as a function of the resonance strength and energy 
under the narrow resonance approximation. From the plot, we see that the 7 Li abundance 
is reduced by 50% for (2 J + l)r e g = 240 keV. This state has J = 5/2 and therefore, a 
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Figure 3: As in Fig. [I] for the resonances in the 9 Be compound nucleus. 
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Table 3: As in Table [I] listing the potential resonances in 9 B. 

value r o fj = 40 keV or more will have substantial impact on the problem. Furthermore, as 
> r c g, we require > 40 keV. This result confirms the conclusion of [16]. Later in 

§ [5] we will see how this compares with theoretical limits. As the decay widths are largely 

unknown, experimental data on the width are needed. 

The state at 17.076 MeV corresponds to a resonant energy of E Tes = 586 keV and is 

beyond the scale shown in Fig. HI A solution using this state is very unlikely. 

4.3 A = 10 Compound Nucleus 

Table HJshows that the 10 Be compound nucleus has energy levels at 17.12 and 17.79 MeV [37] 
which are close to the initial state 7 Li + t at 17.2509 MeV. The former is far below threshold 
and does not contribute to 7 Li destruction. The 17.79 MeV level is around 540 keV above 
the entrance energy and its spin and parity are unknown. The total width [37] is r to t = 112 
keV which implies a small overlap with the entry channel which renders this resonance 
insignificant despite having a number of n exit channels with both ground state and excited 
states of 9 Be. As seen in Fig. [51 the effect of 7 Li + t is small for the interesting region of 
parameter space. 

The 10 B compound nucleus has energy levels at 18.2, 18.43, 18.80 and 19.29 MeV, which 
we investigate. The 18.2 MeV level is uncertain experimentally [38] as indicated in Table EJ 
and hence ideal for parametrizing. There is a 3 He entrance channel a little over 400 keV 
below this level. The current total experimental width is 1.5 MeV which is very large and 
the branching ratios are unknown. The current uncertainty in the E ies is 200 keV. However, 
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Figure 4: As in Fig. [TJfor the resonances in the 9 B compound nucleus. The vertical dashed 
line at 220 keV indicates the experimental central value of the resonance energy of the 16.71 
MeV level. 
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112 ± 35 keV 


7 

n + 9 Be 
n + 9 Be* (1.684 MeV) 
n + 9 Be* (2.4294 MeV) 
n + 9 Be* (2.78 MeV) 
n + Be (3.049 MeV) 
n + 9 Be* (4.704 MeV) 
n + 9 Be* (5.59 MeV) 
n + 9 Be* (6.38 MeV) 
ra + 9 Be* (6.76 MeV) 
n + 9 Be* (7.94 MeV) 
t (elastic) 


3 + 2 eV 

< 77 keV 

< 77 keV 

< 77 keV 

< 77 keV 

< 77 keV 

< 77 keV 

< 77 keV 

< 77 keV 

< 77 keV 

< 77 keV 
78 keV 



Table 4: As in Table [T] listing the potential resonant reactions in 10 Be. 

according to the plot in Fig. [51 even a 200 keV reduction in E res would not be sufficient to 
cause any appreciable destruction of 7 Li as this reaction has negligible effect on the mass-7 
abundance. This is another illustration of the fact that reactions involving direct destruction 
of 7 Li are unimportant. 

The 18.43 MeV level is better understood |39j and with a resonance energy of ~ 640 keV 
for the 7 Li + 3 He initial state (Table [6]) and a total width of 340 keV has a lower entrance 
probability and therefore is likely to be ineffective. This is evident from Fig. |5j This level 
is also a sub-threshold resonance for the 7 Be + t state (Table [7]), with resonance energy, 
E res — —239.1 keV. This is far below threshold rendering it ineffective. 

Staying with 7 Be + 1, the closest energy level above the entrance energy of 18.669 MeV 
is the 18.80 MeV (2 + ) level (Table |8]), which corresponds to a resonance energy of ~ 130 
keV [38J . The exit channel widths for p and 3 He are unknown experimentally and thus, this 
is a candidate for parametrization. There is a weak upper limit on r to t < 600 keV (38] , which 
for J* = 2 is off scale in Figs. [6] and [71 The contour plot in Fig. [6] show that for a central 
value of resonance energy of ~ 130 keV shown by the vertical dashed line, resonance strength 
of just under 1 MeV is required which is very high. Also, parity requirements force L = 1, 
which will cause suppression of this channel. We note that there is no quoted uncertainty for 
this energy level and neighboring levels have typical uncertainties of 100-200 keV. Therefore 
it may be possible (within 1-2<t) that the state lies at an energy of 100 keV lower and would 
energetically, have a chance at solving the 7 Li problem. This is true for the p exit channel. 
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Figure 5: As in Fig. [I] for the resonances in the 10 Be compound nucleus involving the initial 
state 7 Li-K (left), and in the 10 B compound nucleus involving the initial state 7 Li+ 3 He. 
( right) . 



The 3 He exit channel may also reduce mass 7, through the formation of the 7 Li which is 
much easier to destroy. This is reflected in Fig. [71 which shows that at resonance energies of 
< 100 keV, a strength of a few 100 keV, but less than 600 keV may be sufficient to achieve 
comparable destruction of 7 Be as the 16.71 MeV resonance. The caveat is that for such values 
of strengths, the narrow resonance approximation does not hold true and this may lead to a 
reduced effect. Nevertheless, this is yet another case deserving a detailed comparison with 
the theoretical limits which will follow in § [5j Once again, definitive conclusions can be 
drawn only based on experimental data. 

The 19.29 MeV level (Tabled]) is energetically harder to access and with a total width of 
only 190 keV, it is unlikely to be of significance, despite being less studied. 

The 10 C nucleus |40] appearing in Table [10] shows large uncertainties and experimental 
gaps at higher energy levels which may be relevant to entrance channels involving 7 Be. 
Reactions involving the 7 Be + 3 He initial state could contribute in destroying 7 Be if there 
exists a resonance in the parameter space shown in the Fig. [SJ These reactions win over 
those involving the 7 Be + t state, because 3 He is substantially more abundant than t, but 
are worse off due to a higher Coulomb barrier. The entrance energy for 7 Be + 3 He is 15.0 
MeV. As one can see from the figure, a 1~ or 2~ state with a resonance energy of either -10 
keV or 40 keV corresponding to energy levels of 14.99 and 15.04 MeV respectively with a 
strength as high as a few 10 's of keVs is what it will take to solve the lithium problem with 
this initial state. Thus, any 10 C resonance near these energies which may have been missed 
by experiment may be interesting as a solution to the lithium problem; we return to this 
issue in more detail in § [5j 
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Figure 6: As in Fig. [T]for the resonances in the 10 B compound nucleus involving initial states 
7 Be+t. 
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Figure 7: As in Fig. [I] for the resonances in the reaction 7 Be(t, 3 He) 7 Li. 
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Table 5: As in Table IT] listing the ground and excited final states for the 18.2 MeV energy 
level in 10 B. 

4.4 A = 11 Compound Nucleus 

For n B [41J, Table [TUl shows that the entrance channel, 7 Li + a is at 8.6637 MeV which is 
103.7 keV above the resonant energy level at 8.560 MeV and ~ 260 keV below the resonant 
energy level at 8.92 MeV. Parity demands angular momentum to be 0. Both states are at 
relatively large |-E res | and are not capable of making a sizable impact on the 7 Li abundance. 
Table [TUl further lists states 9.19 MeV (which requires L = 3 and has a total width of < 2 eV) 
and 9.271 MeV (whose decay is dominated by the elastic channel) which have progressively 
larger resonant energies and are unlikely to provide a solution. 

For 11 C [12], the entrance channel, 7 Be + a is at 7.543 MeV which is 43 keV above the 
resonant energy level at 8.560 MeV and 557 keV below the resonant energy level at 8.10 
MeV. 

As seen in Fig. [91 we find that the sub-threshold resonance in the 11 C nucleus, produces 
a very insignificant effect on 7 Be in agreement with the claim in [16] . The super-threshold 
resonance states are also too far away at resonance energies, 557 keV and 260 keV for 
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Table 6: As in Table [T] listing the ground and excited final state channels for the 18.43 MeV 
energy level in 10 B for the 7 Li + 3 He initial state. 

7 Be(a,7) n C and 7 Li(a,7) 11 B respectively. 

However, Fig. [9] shows that the presence of a (missed) resonance at resonance energies of 
few tens of keV 's, requires a very meagre strength of the order of tens of meV 's to destroy 
mass 7 substantially. Strengths of this order are typical of electromagnetic channels. It is 
difficult to assess the probability that a U C state at 7.55 MeV has been overlooked. 

5 Reduced List of Candidate Resonances 

Having systematically identified all possible known resonant energy levels which could affect 
BBN, we find most of these levels are ruled out immediately as promising solutions, based on 
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p + 9 Be* (6.38 MeV) 


unknown 








3 






p + 9 Be* (6.76 MeV) 


unknown 















p + 9 Be* (7.94 MeV) 


unknown 








2 






p + 9 Be* (11.283 MeV) 


unknown 















p + 9 Be* (11.81 MeV) 


unknown 








1 






d + 8 Be 


unknown 








1 






d + 8 Be* (3.03 MeV) 


unknown 








1 






d + 8 Be* (11.35 MeV) 


unknown 








1 






a + 6 Li 


unknown 








1 






a + 6 U* ( 2.186 MeV) 


unknown 








1 






a + 6 Li* ( 4.31 MeV) 


unknown 








1 






a + 6 Li* ( 5.37 MeV) 


unknown 








1 






a + 6 Li* ( 5.65 MeV) 


unknown 















t (elastic) 


unknown 



Table 7: As in Table [T] listing the ground and excited final state channels for the 18.43 MeV 
energy level in 10 B for the 7 Be + 1 initial state. 



their measured locations, strengths, and widths. As expected, the existing electromagnetic 
channels are too weak to cause significant depletion of lithium owing to their small widths. 

From amongst the various hadronic channels listed in the tables above, we have seen 
that all channels are unimportant except three which are summarized in table [TTJ The 
7 Be + d channels involving the 16.71 MeV resonance in 9 B, the 7 Be + t channels involving the 
18.80 MeV resonance in 10 B and 7 Be + 3 He channels. These are ones where a more detailed 
theoretical calculation of widths is required to decide whether they may be important or 
not. For each reaction, the Wigner limit, eq. ( !T2|) . to the reduced width 7 2 imposes a bound 
on via eq. (TTTTl . Specifically, the penetration factor, Pi(E,a), must be estimated to see 
if the required strengths according to Figs. HJ El El and M to solve the problem are at all 
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a + 6 Li 


unknown 








2 






a + 6 Li* (2.186 MeV) 


unknown 








2 






a + 6 Li* (3.56 MeV) 


unknown 















a + 6 Li* (4.31 MeV) 


unknown 















a + 6 Li* (5.37 MeV) 


unknown 








2 






a + 6 Li* (5.65 MeV) 


unknown 








1 






t (elastic) 


unknown 



Table 8: As in Table [T] listing the ground and excited final state channels for the 18.80 MeV 
energy level in 10 B for the 7 Be + 1 initial state. 

attainable. The penetration factor is given by 

P ^ = Gl { E, a) \n(E,a) (16) 

where Gl{E,o) and F L (E,a) are Coulomb wavef unctions. 

We note that the Coulomb barrier penetration factor decreases as the energy of the 
projectile and/or the channel radius, a, increases. For a narrow resonance, the relevant 
projectile energy is E w E ies , which is set by nuclear experiments (where available) and 
their uncertainties. The channel radius corresponds to the boundary between the compound 
nucleus in the resonant state and the outgoing / incoming particles. Therefore, the channel 
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3 He 


unknown 
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a + 6 Li 


unknown 








1 






a + 6 Li* (2.186 MeV) 


unknown 








1 






a + 6 Li* (4.31 MeV) 


unknown 








1 






a + 6 Li* (5.37 MeV) 


unknown 








1 






a + 6 Li* (5.65 MeV) 


unknown 















t (elastic) 


unknown 



Table 9: As in Table [T] listing the ground and excited final state channels for the 19.29 MeV 
energy level in 10 B for the 7 Be + 1 initial state. 

radius depends on the properties of the compound state and the particles into which it 
decays. 

Consider the case of 7 Be + <i, which has resonance energy, E res = 220 ±100 keV and initial 
angular momentum, L init = 1. A naive choice for the channel radius is the "hard-sphere" 
approximation, 

a 12 = 1.45 (A 1 / 3 + ^ /3 )fm (17) 

which gives 027 = 4.6 fm. Using the Coulomb functions, T\ is of order a few keVs. The 
corresponding strength, r e g should be essentially the same and we further gain a factor of 
6 from the spin of this state. This suggests by using figure HI that this resonance should fall 
short of the width required to solve or even ameliorate the problem. 

However, reactions involving light nuclides including A = 7 are found to have channel 
radii exceed the hard-sphere approximation [16]. We thus consider larger radii and find that 
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Table 10: As in Table [I] listing resonances in 10 C, n B and n C. 



for values higher than around 10 fm, we get a width which has the potential to change the 
7 Li abundance noticeably. The Wigner limit 



3tf 
2 ii E, 



a 2 = (1{ 



gives a larger radius, 027 = 13.5 fm, which gives one a better chance of solving the problem. 
This is consistent with the conclusions drawn by [16]. 

For the 7 Be + t initial state, the 18.80 MeV state of 10 B has a resonance energy E = 
0.131 MeV and Lj n # = 1. There is no experimental error bar on the resonance energy. The 
hard sphere approximation gives a 37 = 4.9 fm. This gives a width, r 1; which is less than 
a tenth of a keV, and is orders of magnitude lower than what is needed. In the spirit of 
what we did in the earlier case, using eq. (|T8|) gives a channel radius, a 37 = 15 fm improving 
the situation by almost 2 orders of magnitude in l\. If, in addition to increasing 037, the 
resonance energy were to be higher by 100 keV, then l\ could be large enough to change 
the 7 Li abundance noticeably. 

The 7 Be + 3 He initial state will have 10 C as the compound state. The structure of the 
10 C nucleus is not well studied experimentally [10] nor theoretically. In particular, we are 
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Figure 8: As in Fig. [T]for the resonances in 10 C involving initial state 7 Be + 3 He 



unaware of any published data on 10 C states near the 7 Be + 3 He entrance energy, i.e., states 
at or near E ex ( 10 C) ~ Q( 7 Be + 3 He) = 15.003 MeV. To our knowledge, there has not been 
any search for narrow states in this region. The potential exit channels of importance are 
9 B +p and 6 Be + a. Because J^( 3 He) = 1/2+ J^( 7 Be) = 3/2~ to have L init = and thus 
no entrance angular momentum barrier would require the 10 C state to have 

r = (1 or 2)- (19) 

Because J 7r ( 9 B) = 3/2 _ , the entrance channel spin and parity required to give L init = will 
also allow Lg n = for the 9 B+p. On the other hand, in the final state 6 Be + a both 6 Be and 
4 He have J n = + . Thus if the putative 10 C state has J 71 " = l - , this forces the 6 Be + a final 
state to have Lg n = 1, and thus this channel will be suppressed by an angular momentum 
barrier relative to 9 B + p. 
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Figure 9: As in Fig. [T]for the resonances in n C involving initial states 7 Be+a. 
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Using eq. (ITT)) , we again get 037 = 4.9 fm. Taking L init = and E = 0.2 MeV, T is 
about 10~ 3 keV and is extremely small. However, the penetration factor is highly sensitive 
to the channel radius and a relatively small increase in a increases the width by orders of 
magnitude. Increasing the energy does reduce the penetration barrier, but a higher width 
is required due to thermal suppression. In order to get a sizable width, which is required to 
solve the problem according to figure [U a 37 must be > 30 fm. At this energy, this radius is 
somewhat larger than what is afforded by Eq. [T8j 
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unknown 
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7 Be + t 


1 


1 
1 

2 


130.9 keV 


< 600 keV 


p + 9 Be* (11.81 MeV) 
:i He 
a 


unknown 
unknown 
unknown 


10 C, 
unknown 


7 Be + 3 He 


unknown 


unknown 
unknown 


unknown 
(Q = 15.003 MeV) 


unknown 


P 

a 

3 He (elastic) 


unknown 
unknown 
unknown 



Table 11: This table lists surviving candidate resonances. 



6 Discussion and Conclusions 

The lithium problem was foreshadowed before precision cosmic microwave background data, 
was cast in stark light by the first-year WMAP results, and has only worsened since. While 
astrophysical solutions are not ruled out, they are increasingly constrained. Thus, a serious 
and thorough evaluation of all possible nuclear physics aspects of primordial lithium pro- 
duction is urgent in order to determine whether the lithium problem truly points to new 
fundamental physics. 

Reactions involving the primordial production of mass-7, and its lower-mass progeni- 
tor nuclides, are very well studied experimentally and theoretically, and leave no room for 
surprises at the level needed to solve the lithium problem [21[TJ3[T75] . Lithium destruction re- 
actions are less well-determined. While the dominant destruction channels 7 Be(n,p) 7 Li and 
7 Li(jo, ot)a have been extensively studied, in contrast, the subdominant destruction channels 
are less well-constrained. 

We therefore have exhaustively cataloged possible resonant, mass-7 destruction channels. 
As evidenced by the large size of Tables [THTU1 the number of potentially interesting compound 
states is quite large. However, it is evident that the basic conservation laws such as angular 
momentum and parity coupled with the requirement of resonant reactions to be 2-3 times 
the 7 Be(n,p) 7 Li rate prove to be extremely restrictive on the options for a resonant solution 
to the lithium problem, and reduces the possibilities dramatically. 

Given existing nuclear data, there are several choices for experimentally identified nuclear 
resonances which come close to removing the discrepancy between the lithium WMAP+BBN 
predictions and observations as tabulated in § O The 16.71 MeV level in 9 B compound 
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nucleus, and the 18.80 MeV level in the 10 B compound nucleus are two such candidates. It 
is possible, however, that resonant effects have been neglected in reactions passing through 
states which have been entirely missed. In all of the plots above, we have illustrated the 
needed positions and strengths of such states, if they exist. One possibility involving the 
compound state 10 C is poorly studied experimentally, especially at higher energy states close 
to the Q- value for 7 Be + 3 He. 

Any of these resonances (or a combination) could offer a partial or complete solution 
to the lithium problem, but in each case, we find that large channel radii (a > 10 fm) are 
needed in order that the reaction widths are large enough. We confirm the results of Cyburt 
and Pospelov [16] in this regard concerning 7 Be + d, and we also find similar channel radii 
are needed for 7 Be + t, while larger radii are required for 7 Be + 3 He. Obviously, nature 
need not be so kind (or mischievous!) in providing such fortuitious fine-tuning. But given 
the alternative of new physics solutions to the lithium problem, it is important that all 
conventional approaches be exhausted. 

Thus, based on our analysis, quantum mechanics could allow resonant properties that can 
remove or substantially reduce the lithium discrepancy. An experimental effort to measure 
the properties of these resonances, however can conclusively rule out these resonances as 
solutions. If all possible resonances are measured and found to be unimportant for BBN, 
this together with other recent work [15J, will remove any chance of a "nuclear solution" 
to the lithium problem, and substantially increase the possibility of a new physics solution. 
Thus, regardless of the outcome, experimental probes of the states we have highlighted will 
complete the firm empirical foundation of the nuclear physics of BBN and will make a crucial 
contribution to our understanding of the early universe. 

We are pleased to acknowledge useful and stimulating conversations with Robert Wiringa, 
Livius Trache, Shalom Shlomo, Maxim Pospelov, Richard Cyburt, and Robert Charity. The 
work of KAO was supported in part by DOE grant DE-FG02-94ER-40823 at the University 
of Minnesota. 



A The Narrow Resonance Approximation 

Consider a reaction A + b — > C* — > c + D, which passes through an excited state of the 
compound nucleus C*. We treat separately normal and subthreshold reactions, defined 
respectively by a positive and negative sign of the resonance energy E res = E cx — Qc, where 
E ex is the excitation energy of the C* state considered, and Qc = + A(B) — A(C*). 

In general, the thermally averaged rate is 

. , [d 3 v e-» v2 ' 2T av [%~ m , /2 f 00 , „ „ , > FIT . , 

J<Pve-^ I dEE ^- E ' T < 20 > 

For a Breit-Wigner resonance with widths not strongly varying with energy, this becomes 

<"> = (2,^ X d V-£U 2 + (rW2)* (21 > 
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Thus the thermal rate is controlled by the integral of the Lorentzian resonance profile mod- 
ulated with the exponential Boltzmann factor. 

The narrow resonance approximation has usually only been applied to the normal reso- 
nance case, and assumes that the total resonance width is small compared to the temperature: 

r tot < t. 



A.l Narrow Normal Resonances 

In the normal or "superthreshold" case, the integral includes the peak of the Lorentzian 
where E = E ves . The narrow condition then guarantees that over the Lorentzian width, the 
Boltzmann factor does not change appreciably, and so we make the approximation 

exp(-|) «exp[-|) (22) 



TJ IT 

where we choose the "typical" energy to be the peak of the Lorentzian, E = E res . Then the 
integral becomes 



Furthermore, it is usually also implicitly assumed that the resonance energy is large compared 
to the width: E Tes 3> r tot . Then the integral gives 27r/r tot , and the thermally averaged cross- 
section under this approximation is given by 



2 \ 3/2 

Mr tot «T, Srcs = ^r cff (— e~ E -' T (24) 



jxT J 

= 2.65 x l(r 13 /i- 3/2 to T eS T~ 3/2 exp(-l 1.605 E ICS /T 9 ) cmV 1 (25) 

where the latter expression has T 9 = T/10 9 K. 

Note however, that eq. (|2~T!) is exactly integrable as it stands and does not require we 
make the usual E res ^> r tot approximation. Thus for the normal case we modify the usual 
reaction rate and instead adopt the form 

\ UVJ narrow, normal (""^')rtot < S7 1 ,i? res /(2£ res /r tot ) (26) 

Here we introduce a temperature-independent correction for finite E rcs /Y tot (still with E res > 
0) 

f(u) = — I — arctanw . (27) 

2 7T 

This factor spans / — > 1/2 for E res T tot to / — > 1 for E Tes ^> r tot . 

In practice, we adopt a slightly modified version of the correction factor in our plots. 
Recall that in Figs. I2H9J we show results for lithium abundances in the presence of resonant 
reactions with fixed input channels, but without reference to a specific final state. Without 
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the correction factor, the resonant reaction rate is characterized by two parameters, E ves and 
r eff . These two parameters are insufficient to specify the correction factor, which depends 
on E TCS /r to t- Rather than separately introduce r tot , we instead approximate the correction 
factor as f(2E TCS /T e s). Because T c s < r tot and / is monotonically increasing, this always 
underestimates the value of / and thus conservatively understates the importance of the 
resonance we seek (but the approximation is never off by more than a factor of 2 in the 
normal case). 



A. 2 Narrow Subhreshold Resonances 

Still making the narrow resonance approximation r to t "C T, we now turn to the subthreshold 
case, in which E res < 0. To make effect of the sign change explicit, we rewrite eq. (12 ip as 

, , _ CoT init r fin f°° e~ E/T 

{(TV) ~ 2(2^)3/2 J dJ "(E + \E ies \Y + (T tot /2Y 

Now the integrand always excludes the resonant peak, and only includes the high-energy 
wing. As with the normal case, the narrowness of the resonance implies that the Boltzmann 
exponential does not change much where the Lorentzian has a significant contribution, and so 
we again will approximate e~ E l T ~ e~ E l T . Since we avoid the resonant peak, the choice of E 
not as straightforward in the subthreshold case where we took E = E rcs . This choice makes 
no sense in the subthreshold case, because the e~ Elcs ^ T > 1 in the subthreshold case, yet 
obviously kinetic energy E > and thus the Boltzmann factor must always be a suppression 
and not an enhancement! 

Yet clearly \E res \ remains an important scale. Thus we put E = u\E TeB \, and we have 
examined results for different values of the dimensionless parameter u. We find good agree- 
ment with numerical results when we adopt u ~ 1, i.e., E = \E Ie8 \. Thus for the subthreshold 
case we adopt a reaction rate which is in closely analogous to the normal case: 

/ 2tx \ 3 ^ 2 

(c^narrow, subthreshold = W T eS ( — ^ J e ' Efes ^ T f (— 2 1 E res \ /T tot ) (29) 

Similarly to the normal the reaction becomes increasingly off- resonance, i.e., as 

\E Tes \ grows, there is an exponential suppression. In addition, the correction factor has 
limits / — ¥ 1/2 for |-E res | "C r tot , and / — ¥ as |-E res | S> r tot . Finally, note that, as a function 
of .E^es, our subthreshold and normal rates match at E res = 0, as they must physically. 
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